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Abstract. In this paper we investigate the properties of Bose gases with 
Raman-induced spin-orbit (SO) coupling. It is found that the SO coupling can 
greatly modify the single particle dcnsity-of-state, and thus lead to non-monotonic 
behavior of the condensate depletion, the Lee-Huang- Yang correction of ground- 
state energy and the transition temperature of a non-interacting Bose-Einstein 
condensate. The presence of the SO coupling also breaks the Galileaan invariance, 
and this gives two different critical velocities, corresponding to the movement of 
the condensate and the impurity respectively. Finally, we show that with SO 
coupling, the interactions modify the BEC transition temperature even at Haxtree- 
Fock level, in contrast to the ordinary Bose gas without SO coupling. All results 
presented here can be directly verified in the current cold atom experiments using 
Raman laser-induced gauge field. 
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1. Introduction 

Recently one of the major progresses in cold atom physics is the realization of spin- 
orbit (SO) coupling for neutral atoms with the idea of light induced synthetic gauge 
fieldpQ [21 [21 13]. Previously, effects of spin-orbit coupling have only been studied 
in fermionic matters such as electron gases, while this progress opens up the new 
opportunity of studying SO coupling in bosonic systems. This investment is interesting 
because SO coupling changes many basic properties of Bose-Einstein condensation 
(BEG) in both quantitative and qualitative way. 

Although SO coupling emerges at single particle level, it may have a significant 
effect to the many-body behavior due to the modification of the single particle 
density-of-state (DoS). One particular example is the Rashba SO coupled Bose 
gases. Because of the dramatic change of low energy DoS, many intriguing many- 
body phenomena have been predicted there, and some questions still have not been 
thoroughly understood yet. While a number of theoretical works have studied its 
low-temperature properties, including ground state phases, quantum and thermal 
fluctuations, Bose condensation transition and superfluiditiy [3 [3 [H [H [TOl [TT] . 
the experimental realization of Rashba SO coupling remains to be a great challenge. 
On the other hand, the type of equal Rashba and Dresselhaus SO coupling has been 
realized in current experiment with two photon Raman transition. Although the 
change of low-energy DoS is not as dramatic as Rashba case, such a SO coupling leads 
to a non-monotonic behavior of low-energy DoS as a function of Raman strength, which 
should result in nontrivial manifestation of many-body properties. So far, a few papers 
have studied the ground state phase diagram, collective modes and response function 
for such SO coupled condensate [T^l [T31 [T31 [ISl UHl [IZ]: but there are still several 
fundamental properties which have not been explored, in particular, the depletion of 
the condensate, the BEG transition temperature and the superfluid critical velocities. 
In this work we report our theoretical studies of these properties. 

2. Single particle Hamiltonian and mean-field phase diagram 

2.1. Single particle Hamiltonian 

In current experiments, two counter-propagating Raman beams are applied to couple 
two hyperfine levels of an alkali atom [H l2l [3l|4] , which gives rise to the single-particle 
Hamiltonian as [THl (set h = 1) 

^ _( (p-Ke^f /2m + S/2 Q./2 

where kr is the recoil momentum, U, is the Raman coupling strength, and 5 is the two- 
photon detunning. With Pauli matrix, the Hamiltonian in Eq. ([!]) can be rewritten 
as 

p2 

Hq = Bp ■ (T + Er (2) 

im 

where Er — /c^/(2m) is the recoil energy, and Bp = {—^/2,0, krPx/rri — d/2) 
depends on momentum and yields the locking between spin and momentum. This 
Hamiltonian preserves spatial translational symmetry, and momentum p is a good 
quantum number. Another quantum number of this Hamiltonian is "lielicity" /i = ±, 
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Figure 1. The phase diagrams of *^Rb (a) and ^''Na (b). The two pseudo- 
spin states are m = 0,-1 states of F = 1 hyperfine states. On the right of 
red dashed line, the single particle spectrum has one local minimum, and there 
is only one plane wave phase (PW3). On the left of red dashed line, the single 
particle spectrum has two local minima, and the horizontal blue line separates 
two plane- wave condensates at different momenta (PWl and PW2). Within the 
region rounded by two blues lines it is the stripe phase, in which the condensate 
coherently occupies two different momenta. Here Er = 2n X 2.2kHz. 



which denotes for Bp parallel or anti-parallel to spin cr. Thus, the eigen-energies of 
two helicity branches are given by 
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and their wave-functions are given by 
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^JiPxkr/m-5/2f + ^1^4^ 
With the spectrum, the single particle DoS can be calculated straightforwardly as 

^(^) = t7 E - ^P^+) + '^d(£ - ep._)] (5) 



where (5d is Dirac delta function. 



2.2. Mean-field phase diagram 

The interaction of a two-component Bose gas is generally written as 

j3 



1 
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With mean-field approximation, the interaction energy is given by 

^i = \j ^'"^ (5n"r(^) + aunliv) + 2g^^n^{r)n^{r)) (7) 

When 5 and Vl are both small, e_ has two local minima located at k±. Without loss 
of generality, we can assume the condensate wave function as 
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where < a < 7r/2, cos a and sin a are superposition coefficients. It gives 

— no [cos^ a cos^ 0+ + sin^ a cos^ 6^ + sin 2a cos 6*+ cos 6'_ cos((5p • r)] , 

ni — no [cos^ a sin^ + sin^ a sin^ 0_ + sin 2a sin 6*+ sin 6^ cos((5p • r 

where dp = p+ — P-- In practices, one can straightforwardly substitute the expression 
of density Eq. [9] into the interaction energy Eq. [Tj and the interaction energy £j 
becomes a function of a, p± and 9±. Then by minimizing the total energy respect to 
a, p± and 9±, one can finally obtain the ground state condensate wave function for 
given fi, 6 and interaction parameters g^^, and g^^. Thus, a phase diagram can 
be constructed. More details have been discussed in Ref. [TJ [T^l US] and here we will 
just emphasize some general features. 

(1) If the energy minimization gives a = or a = 7r/2, there is only one 
momentum component in the condensate wave function, and the density for both 
spin component are uniform, as it can be easily seen by setting a = or 7r/2 in Eq. 
([9]). This is the "plane wave" condensate. While if the energy minimization gives 

< a < 7r/2, both ^(r) and %(r) have spatially periodic modulation. This is named 
as "stripe" condensate [5 . 

(2) Because of SO couphng, the densities of both components depend on 
momentum. If the interaction is SU{2) invariant, i.e. g^-f = g^ = g-f^ = g, the 
interaction energy becomes £1 = {g/2) J d'^r{n^{r) + n^(r))^. In this case, if the 
condensate is in the "plane wave" phase, the interaction energy is independent of 
momentum, and one can take p± as k± which minimizes the single particle energy. 
However, because the two pseudo-spin states are taken as two hyperfine level of atoms, 
they do not have to obey SU{2) spin rotational symmetry. Thus, for a general case, 
three interaction parameters g-f^, g^^ and g^j^ are all unequal. Therefore, even for the 
plane wave condensate, the interaction energy also depends on momentum p±. That 
is equivalent to say, even at mean-field level the self-energy correction has momentum 
dependence, which effectively modifies the single particle dispersion and changes the 
location of its minimum. In contrast, without SO coupling, the mean-field self-energy 
correction is just a constant shift of single particle energy. Thus, as first pointed out 
by Ref. [T3] , p± is shifted away from single particle minimum k± [|] 

(3) The stripe phase has distinct properties from plane-wave phase. For the 
stripe phase, the total density will have spatial modulation, because the spin wave 
function at p+ and p_ are not orthogonal, i.e., cos^p^ cos6'p_ +sin6'p^ sin6'p_ 7^ Cj§] 
Thus, in this system the stripe phase is not favored by density interaction part 
{g/2) J (Pr{n^{r) + n^(r))^. In another word, if the interaction is SU{2) invariant, 

1 However, for some special case, such as pure Rashba SO coupling without Zeeman field considered 
in Ref. [5|, the mean- field self-energy correlation is also momentum independent. 
§ This is also different from pure Rashba case where the spin wave function are orthogonal for 
opposite momentum, and the total density is uniform for stripe phase 
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Figure 2. Upper panel: Single particle dispersion £±{px) as a function olpx/kr; 
Lower panel: DoS as a function of e— Emin for f2 = ^Er, AEr and &Er, respectively. 
Epmin denotes the minimum value of single particle energy. The dashed lines are 
DoS for the case without Raman coupling. 

stripe phase will not exist in this system. The difference in 5-^^, g^^ and are 
necessary for stabilizing the stripe phase. Moreover, since the non-uniform term in 
total density increases as increases. Thus, the stripe phase, if exists, should be 
found in small f2 regime of the phase diagram. 

(4) Consider the limit f2 = 0, if g^^g^^ — <?^^ > 0, a homogeneous mixture of two 
components is stable against local density fluctuations, and there will be a mixed phase 
within certain detuning window. Such a mixed phase will turn into stripe phase once 
Q, becomes non-zero, for instance, for *^Rb case in Fig. 1(a). While if — .g|^ < 0, 

the mixed phase is not stable against phase separation even for zero 17, and there will 
be no stripe phase in the phase diagram, for instance, for ^^Na case in Fig. 1(b). 

Hereafter, we should focus only on SU (2) invariant interaction. This is relevant 
for experiment with Rb or Na, because the difference in g^^^, and g^^ are smaller 
than 1%. The generalization to non-S'C/(2) interaction is straightforward. Besides, we 
only consider the plane wave phase, because in these systems, the stripe phase either 
occupies a very small regime of phase diagram or does not exist. 

2.3. Zero-detuning case 

In this work we will particularly focus on the case with 5 — Q for following two reasons. 
(1) Density- of- State Effect. When f2 < 4i?r, t-{kx) has two minima at fc± = 

-iikr\j 1 — ^TE~^ ^ ^^"^ when > 4£'r, e-{kx) has one single minimum at — 0. 
Expanding the dispersion around its minimum, one gets the effective mass in x 
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direction as 

m* = l -(l-ii?)" ^<4^'- . (10) 
[ m (1 - ^) n > AEr 

Hence, the low-energy DoS increases with fl when < AEr and decreases with 51 
when rt > AEr, as shown in Fig. [2j The most intriguing point is at £7 = AEr when the 
single particle dispersion behaves as ~ p!^ at the lowest order and the low-energy DoS 
reaches its maximum. As we shall see in later discussion, this has important physical 
consequences in the superfluid critical velocity and the BEC transition temperature. 

(2) Z2 Symmetry and Magnetization. When 57 < AEr, bosons condense into one 
of the minima, which breaks the Z2 symmetry. The Bose condensate will have finite 
magnetization. While when D, > AE^., bosons condense at zero- momentum state and 
the condensate is non-magnetic. Thus, there will be a magnetic phase transition at 
f2 = AEr associated with the Z2 symmetry breaking, and a divergent spin susceptibility 
has been predicted and experimentally found [H [15] . We note that such a transition 
exists only for 5 = 0, since for non-zero S the Hamiltonian does not possess the Z2 
symmetry, and the condensate phase is always magnetic. 



3. Bogoliubov Theory and Superfluid Critical Velocity 

3.1. Bogoliubov Spectrum 

We study the fluctuations around the condensate with Bogoliubov theory. Considering 
a plane- wave condensate at momentum po = poe^, the field operator can be expanded 
as 

V'(r) = ^(r) + (5i/'(r) (11) 
where <p(r) is the condensate wave- function 

^(r) = ( ) exp i^pox) (12) 

and satisfies the Gross-Pitaevskii (GP) equation 

Ho(»)H-.H(::i.)^.( ™:?.)- w 

Defining #t = (■0p„+q^^, ^p^+q^, V'po-q,T' ^po-q,;)' the Bogoliubov Hamiltonian 
for the fluctuation part can be written as 

^ - I] *q^q*q -^Yl [^po-q.t + ^Po-q,; - 2^ + ^gm] 

where 

""^K Xo(po-q) + S^ ) ^'^^ 

Kq ~ Hq — /I is the grand Hamiltonian of non-interacting system, Sn and Ea are 
normal and anomalous self-energy respectively 

V - nr, f sin^ ^'po + 2 cos2 6lp„ -sin6lpocos6'p„ X (,.^ 
-S'^oy _ ^^g2 ^ 2 sin2 Op, j ' 
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Figure 3. Bogoliubov spectrum for Q/Er = 2,4 and 6 for (a-c). gn/Er is fixed 
to be 0.5. 
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Hence the Bogoliubov spectrum is determined by 
" Koipo + q) + Y.^ ~ E 



(16) 



Det 



(17) 



-Sa -ifo(po-q)-SN-S 

Because of the emergence of off-diagonal long range order, the excitation energy 
should be gapless in the long wave-length limit q — ^ 0, which requires 



Det 



Ya{Kq{po) + Sn) 'Sa - (Koipo) + Sn) 







(18) 



This equation is indeed satisfied because GP equation for the condensate wave-function 
can be rewritten as 



Det 



Sa 



= 



(19) 



Eq. (191 can also be regarded as a generalization of Hugenholtz-Pines relation 
/Li = Sn — Sa to the SO coupled Bose gases [19 . 



Solving Eq. (17) gives rise to the entire Bogoliubov spectrum i?q,±. Examples 
are displayed in Fig. 3] for various Q/Ei-. For q^^ —> 0, the low-energy excitation is the 
phonon mode with a linear dispersion E_{qx) = sq^- Along the direction of Raman 
beam, the sound velocity is given by 



w rn V 



1 



f]2 



16E? J 



(20) 



when < AE^ and 



AE, 



(21) 



'^(1 _ 
m \ il 

when 17 > AE,-. These results are consistent with the effective mass approximation 
that gives s = ^ gru^jra* . At = 4£'r, the effective mass diverges, and the low- 
energy phonon mode is quadratic in q^, as shown in Fig. [sj^b). This is a significant 
difference compared to the conventional phonon mode, which is always linear in q. 
When O < 4i?i., the Bogoliubov spectrum has another local minimum at finite q^i 
which is due to the double degeneracy in single particle spectrum. This is attributed 
as "roton minimum" by Ref. |15) . 
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Figure 4. (a) Total depletion as a function of fl/Er for T = 0, O.OSTc and 
O.lTc, where T,, is the condensation temperature (see Section 4). (b) Thermal 
depletion as a function of T/Tc for various values of Q/Er (take gn/Er = 0.5 and 
n/k^ = 1). 



3.2. Quantum and Thermal Condensate Depletion 

The Bogoliubov Hamiltonian can be straightforwardly written as 

+ \Y1 [^-q^+ + -^-q.- " ^Po-q,t - ep„_q,; + 2n- 3gno (22) 

it ~t ',t _ ; 

here ^I'q = (V'q,+ 7 V'q,-: V'-q.+i ^_q,_) = -'V^q "fq are quasi-particle operators, and Mq 
is the Bogohubov transformation matrix which is obtained from eigen equation (17). 
The quantum and thermal depletion of condensate ncx,cr = y X^q^^o^^po+q-o-^Pn+q,"') 
thus can be readily calculated using Mq and £'q,±- 

The depletion fraction as a function of Raman coupling strength D, and 
temperature T is shown in Fig. |4] The contribution to depletion is mainly from 
two parts: the low-energy phonon part and the roton part. At zero-temperature, the 
contribution is dominated by the phonon part because the roton part has a finite 
excitation gap. Thus, due to the non-monotonic behavior of sound velocity discussed 
above, when fl < AEr, the quantum depletion Uex/n increases with J7, reaches a 
maximum ai Q — 4:Er, and then decreases as fl increases. At finite temperature, 
because the roton gap is quite small at small Jl, the roton part will give an significant 
contribution. For < AEr, the contribution from phonon part always increases with 
fl because the decrease of phonon velocity, while the contribution from roton part 
decreases with because the roton gap increases with il. Due to the interplay of 
these two contributions, for small f2, n^x/n first decreases as fi increase, in contrast 
to zero-temperature case. Then n^^jn increases again with and the peak of r^ox/" 
at = A,E^ retains. 

Another interesting manifestation of roton minimum is through the magnetization 
of non-condensed part defined as (jio^^ — 7icxi)/'T'cx- Assuming the condensate 
momentum > 0, the condensate has a positive magnetization (no^ — no^)/no > 0. 
As shown in Fig. [s] (a), at zero-temperature, the depletion also has a positive 
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magnetization (ricxt ~ "■cx^)/?icx > 0. This is because the phonon contribution, which 
dominates the quantum depletion, originates from states whose momenta are nearby 
condensate momentum and possess same magnetization as the condensate. However, 
at finite temperature, for smaU f2, the depletion can have opposite magnetizations 
as the condensate part, as shown in Fig. [sja). And from Fig. [sf^b) one can also 
see that, for instance ai fl — O.bEr, the magnetization soon decreases to negative 
as temperature increases. This is because the roton contribution dominates in this 
regime. Recall that when fl < AEr there is double minimum in this single particle 
spectrum with opposite momentum and magnetization. When condensate takes place 
in one of the minimum, the roton minimum in the excitation appears nearby the 
other minimum. Thus, the contribution from roton minimum displays opposite 
magnetization with condensate part. As we know, the phonon dispersion is linear, 
while the roton dispersion is quadratic. That implies the DoS of the quasi-particle 
near the roton is larger than the one near the phonon part. At finite temperatures 
higher than the roton gap, the thermal depletion near the roton minimum will exceed 
the depletion from phonon modes. This leads to an opposite magnetization of the 
total thermal depletion compared with the condensate. In the further experiment, 
measuring the magnetization of the thermal depletion will thus provide an evidence 
for the exsistence of the roton minimum in the excitation spectrum. 



3.3. Beyond-mean- field correction of ground state energy 

Bogoliubov theory also gives rise to a correction to mean-field energy usually named as 
Lee-Huang- Yang correction [20]. The Lee-Huang- Yang correction can be considered 
as a sum of all zero-point energies of excitation modes at different momenta. Without 
SO coupling, it is a function of n}^^as only and the ratio between energy correction 
to mean-field energy is 

^hBY ^ / 1 N 128 



^ =rLHY(ni/'as) = T^y^. (23) 
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Figure 6. Flhy as a function oi Q/Er for different n^^^as (fixed n/kf = 1). 



While with SO coupling, using Eq. (22 1 one finds that £lhy is given by 



LHY 



2 ^ 



q#0 
H y 

^ fc2 
k 



where the last term is obtained from the renormalization relation of coupling constant 
g via second-order Bohn approximation 
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Therefore Flhy is a function of n^/'^Os, ^/E^ and n/k^, and the results are shown in 
Fig. [6j We find that the Lee-Huang- Yang correction shows a nonmontonic behavior 
as a function of Raman coupling strength and displays a minimum at = A^Ej-. This 
is consistent with the softening of phonon mode at Vic discussed above. The softer 
the phonon mode is, the smaller the zero-point energy will be, which results in the 
minimum in Lee-Huang- Yang correction. 



3.^. Superfluid Critical Velocity 

There are two ways to measure critical velocity [11 . First, the condensate is at rest 
and an impurity moves with finite velocity in the condensate. In this way, one can 
obtain a critical dragging velocity 



Vdrag = min 




where Eq is the quasi-particle energy for momentum q. Critical dragging velocity 
is shown in Fig. 6(a). One can see that vj^ag vanishes at AEr, because the phonon 
spectrum becomes quadratic in at this point and the phonon velocity vanishes. One 
can also see that for fl < AEr, Vdrag is not symmetric for moving along x or along —x. 
This is due to the non-symmetric structure of Bogoliubov spectrum as shown in Fig. 
|3][a). Moving along x, Vdrag is determined by phonon part. Moving along —x, Vdrag 
is determined by the the roton minimum for small VI. Since the roton gap increases 
with VI, the Vdrag along —x also increase with VI untill the roton vanishes. After that 
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Figure 7. Critical dragging(a) and fiowing(b) velocity as a function of Q/Er- 
Here gn = 0.5Er, and Vr = kr/m. 



the Vdrag along —x is also determined by the phonon mode, and will decrease with n 
until fl = AEr. 

Second, one can consider a static impurity and let the condensate move with a 
finite velocity. This determines another critical velocity which can be called critical 
flowing velocity Vflow If the system is Galilean invariant, Vdrag and Vfjow should 
be identical. However, SO coupling breaks the Galilean invariance, and thus these 
two velocities become unequal, as first pointed out in Ref. [11] for Rashba type SO 
coupling. In our system, Galilean invariance is also broken. Considering a Galilean 
transformation in x direction, the Hamiltonian in the moving frame becomes 

Ho (v) = ^ - krCTzf^ + ^i^CFx - v^Px (27) 

= ^ [Pi^ - (krCTz + mvx)f + ^ricTr, - v^krCTz - ]^mv^ (28) 

With a gauge transformation U (x, t) — exp [^imVxX — i (mw^/2) t] , the Hamiltonian 
becomes 

1 2 1 

Ho (v) = ^{Px- Kf^z) + 7,^(^x - V^krCJz (29) 

2m 2 

Compared to the Hamiltonian in the stationary frame, there is an additional velocity 
dependent Zeeman term Vxkr<7z- This term can not be gauged away, and that implies 
the broken of Galilean invariance for a SO coupled particle. The physical effect of this 
term has already been observed in Ref. |4| in collective dipole oscillation experiment. 
Vflow 7^ Vdi-ag is another manifestion of the absence of Galilean invariance. 

To determine Vflow we shall first find out the ground state wave function for the 
Hamiltonian in the comoving frame [Eq. (28l], say 

V9'(x) = ^ ^ - sLV ) '^^^ ■ ^^^'^ 

Then, following similar procedure discussed above, one can find out the Bogoliubov 
spectrum above this new ground state in the comoving frame, denoted by E'^ (v) . 
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Figure 8. Tc (a) and its relative siiift from non-interacting Bose gases ATc/Tc 
(b) of a uniform system as a function of Q/Ei for various interaction strengths. 
Here n/kf = 5. 



Then, observing in the laboratory frame, the excitation spectrum is given by (v) — 
E'^ (v) + V • q. When v is above certain critical value, Eq will start to have negative 
part, which indicates instability of the condensate. This critical values defines Vfjow 
Vflow of this system is shown in Fig. 7(b). We find that, strongly in contrast to Vdragi 
Vflow remains finite across 51 = 4i?r. This is because in the comoving frame with finite 
velocity, the single particle spectrum is always quadratic around its minimum due to 
the velocity dependent Zeeman term, which prevents the softening of phonon mode. 



4. Hatree-Fock theory of normal state and BEC transition temperature 
4.I. Transition temperature of a noninteracting gas 

Now we study the normal state of the system and determine the transition temperature 
of condensate. First, we consider a noninteracting gas. The condensation transition 
temperature of a uniform system is given by (set — 1) 

^%(e-;.)/T_i' (31) 

where chemical potential reaches the bottom of the single particle spectrum, = emin- 
Due to the DoS effect discussed in Section 1, Tc decreases with Q, in the regime 
r2 < 4i?r, reaches a minimum of finite value around 51 = 4£'r, and then increase in the 
regime SI > 4£'r- This non-monotonic behavior is shown in Fig. [8](a). In 51 — >■ limit 
and SI oo limit, one can find a simple relation of transition temperature, 
Tc (SI ^ oo) 

because the low-energy DoS for the later case shrinks to only half of the first one. 
In a harmonic trap, the semiclassical energy of single boson can be expressed as 

£p,± (r) = ^ (P' + fcr) ± ^{krPjmf + -f Knuj^r\ (33) 
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with uj the trap frequency. With semiclassical approximation, the DoS should be 
modified as 

Arap(£) = ^T.J d'r[5,,{e - Sp^r)) + 5,,{e ^ ep,-(r))] (34) 
p 

and the transition temperature can be obtained from 

^-r^ %(f:r;T!'\ (35) 

when chemical potential fi reaches the minimum of single-particle energy at trap 
center. The result is shown in Fig. |9ja). 

In contrast to uniform case, one finds the minimum location of the Tc is shifted 
into smaller ft regime even for a small particle number. While if the particle number 
is large enough, Tc will always monotonically increases with fl. This is because the 
effective mass increases with fi, in a semiclassical sense it will leads an increasing of 
centre density in the trap of the thermal gas[12], ut (0) cx (m*r)'^^^. From a quantum 
view of point, the effective harmonic length a?„ = , of the SO coupled Boson 

no y'rn^muj 

decrease with 51. That also indicates the increasing of the central density. The density 
effect competes with DoS effects, and for large N the former dominates over the latter. 
In 51 — > oo limit, due to the fact that DoS shrinks to one half of that in fl ^ limit, 
the transition temperature in these two regimes follows, 

4-2. Mean-field Shift of Transition Temperature for Uniform Gases 

In the absence of SO coupling, it is well known that the contact interaction between 
the particles does not affect Tc at mean-field level [21], this is because the Hatree- 
Fock self-energy only provides a constant shift of chemical potential. While with SO 
coupling, as shown in the following, the interactions do have a non-trivial effect even 
at mean-field level. 

To construct a self-consistent Hartree-Fock theory, besides the average density 
for each spin component 

"-^ = ^ H ('/'p.^^P.'^ 
p 

we also need to introduce mean-field parameter associating with the spin-fiip term 
1 / J,t 



p 

The spin-fiip term, i^, is due to the Raman coupling between different spin states, which 
breaks the conservation of spin magnetization. In contrast, such a term is absent for 
an ordinary spinor Bose gas without Raman coupling. 

Up to a constant, the mean-field Hamiltonian in Hartree-Fock approximation is 
given by 

p 
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Figure 9. (a) Tc of *^Rb and non-interacting Bose gases in tiie iiarmonic 
trap, (b) Shift of Tc of **^Rb due to interaction as a function of Q,/Er. Here 
Er = 2iT X 2.2kHz. The particle number is Af = 2.5 X 10^, and the trapping 
frequency is tj = 27r X 50Hz. 



Since we are treating the normal state above Tc, spin are always unpolarized with 
= ni = n/2. The mean field Hamiltonian in (371 has the same structure 
as the single particle Hamiltonian, except that the Raman coupling is modified as 
r^cff = + 2gf and the energy zero-point is shifted by a constant. In this sense, we 
can define two dressed helicity branches, with dispersions given by 



1/2 ■ 3 



£p,± = ^ IP' + ^l) + \9n ± ^{krp^/mf + (38) 
Here the mean-field parameter ^ should be solved self-consistcntly in combination with 

^ = ^sin^pcos^p (np^+ - Up _) (39) 
p 

where 9^ is given by Eq. (12. Ih with the bare Raman coupling replaced by rJcff j a-nd 



rip^i are Bose distribution functions with the dispersion £-p,± given by Eq. (38 1. The 
transition temperature Tc is then determined by interaction modified dispersions. Due 
to the non-trivial momentum-dependence of dispersions ( 38 1 , the interactions induce 
a shift of Tc from non-interacting Bose gas. 

Since the occupation of lower helicity branch np._ is always larger than n^.+, f is 
always negative and hence the effective Raman coupling strength ilcft is decreased to 
a smaller value. Based on the non-monotonic behavior of non-interacting Tc discussed 
previously, the minimum of the condensation temperature is shifted by the interaction 
to a larger £7, as shown in Fig. |8] (a). Near the minimum, the correction of Tc change 
its sign rapidly within a narrow region of Q/Er- From Fig. [Sjb), one can clearly 
see that the transition temperature is enhanced for n < AE,- and is suppressed for 
r2 > AE,-, and in between the correction ATc = Tc — Tc shows a most profound effect 
around the Tc minimum. 
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4-3. Interaction Shift of Transition Temperature for Trapped Gases 

Finally we come to discuss the most realistic case where both interaction and trap 
effect are taken into account. We use the Hatree-Fock mean- field theory to include 
interaction effect and use the semi-classical approximation to include trap effect. 

Within the Hatree-Fock and semi-classical approximation, T^ is determined when 
/Lt(T) reaches the minimum of single particle spectrum, and /i(r) is determined by 
particle number conservation 



N 



J d^rn{r,T,^) (40) 



Here the local density n (r, T, fx) is given by 



in which 



and 



£±,k (r) = (k2 -I- kl) ± J{krk^/mf + (r) /4 + V,s (r) , (42) 
2m ^ ' V 



V,s{r) = V{r) + i3/2)gn{r). (43) 

Similar as above, fioff is the renormalized coupling as 

n,s (r) = r! + 2.g^ (r) . (44) 

which also needs to be determined self-consistently as Eq 39 

The numerical solution gives the interaction shift of Tc inside a harmonic trap, 
as shown in Fig. |9] Here we notice that the interaction effect always gives a negative 
shift of ATc, in contrast to the uniform case where ATc can be either positive or 
negative. This is because the presence of the repulsive interaction reduces the central 
density in the trap, leading to a reduced Tc. This effect dominates over the shift 
of r^cff for sufficient large particle number. In Fig. [o] we also plot |ATc|/T° as a 
function of flc, and find the relative shift of reaches a maximum around fl — AE^. 
This can be qualitatively understood from the effective mass approximation. Using 
m*, one can define an effective scattering length a* = as{m* /m) and harmonic 
length aj*Q — aY,o{rn/m*y^^; therefore the shift of Tc in the harmonic trap 23J is 
ATc/T" = -1.32-2^Arl/6 (m* jmfl^. As we have shown in Eqfiol m*/m is 

ho 

maximally enhanced at 17 = 4i?i., so the interaction has the maximum effect to the 
shift of Tc at this point. 



5. Conclusion 



In this paper, we investigate the properties of Bose gases with Raman-induced spin- 
orbit coupling. Our main results are summarized as follows. 

(l)The presence of the SO coupling modifies the single particle spectrum and thus 
the single particle DoS. At = 4i?]., the low energy DoS reaches a maximum and the 
effective mass is maximally enhanced. The direct consequences include the vanishment 
of sound velocity at SI = 4i?i., and the non-monotonic behavior of condensate depletion, 
Lee-Huang- Yang correction of ground-state energy, and the transition temperature of 
a non-interacting Bose-Einstein condensate. 
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(2) The presence of the SO couphng breaks the GaUlean invariance. As a result, 
the critical dragging and flowing velocity, respectively defined in the rest frame of the 
condensate and the impurity, are no longer identical. 

(3) In the presence of the SO coupling, a roton minimum will appear in the 
excitation spectrum in the regime of < 4Er. As a result, the thermal depletion of 
the condensate can possess an opposite magnetization with the quantum depletion. 
Moreover, the critical dragging velocity exhibits asymmetry along different directions. 

(4) In the presence of the SO coupling, the interactions shift BEC transition 
temperature Tc even at a Hartree-Fock level. In both homogeneous and trapped 
systems, the interaction shift of Tc shows maximum around Q — -iE^, where the 
interactions take the largest effect due to the enhancement of DoS and the effective 
mass. 

In conclusion, we have shown that the Bose gases with Raman-induced SO 
coupling can exhibit a number of non-trivial properties, as summarized above. The 
results revealed here can be directly verified in the current cold atom experiments 
using laser-induced gauge field. 
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